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Abstract 

The migration of radionuclides from a borehole repository located about 20 km from the Akwapim fault line which 
lies in an area of high seismicity was analyzed for some selected radionuclides. In the event of a seismic activity, 
fractures and faults could be rejuvenated or initiated resulting in container failure leading to the release of radionuclides. 
A numerical model was solved using a two-dimensional finite element code (Comsol Multiphasics) by taking into account 
the effect of heterogeneities. Results showed that, the fractured medium created preferential pathways indicating that, 
fault zones generated potential paths for released radionuclides from a radioactive waste repository. The results 
obtained showed that variations in hydraulic conductivity as a result of the heterogeneity considered within the 
domain significantly affected the direction of flow. 
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Introduction 

The site for the Borehole Disposal Facility is within the 
Accra Plains and is situated on the boundary between 
the Dahomeyan System and the Togo Series (Figure 1), 
both of which are of Precambrian age. The Dahomeyan 
is the major bedrock formation underlying the site and 
consists of gneiss and schist while the Togo Series are 
predominantly composed of quartzite and phyllite. The 
Akwapim hills and the Togo range mark the line of a 
major active fault zone that runs North-East and South- 
west to the coastal boundary fault (Muff & Efa 2006). 

The whole of the site is covered by loose unconsoli- 
dated and weathered material that may reflect the pres- 
ence of troughs formed by down faulted blocks which 
indicates the existence of seismic activity in the geologic 
past and it probably results from movements along the 
Akwapim fault line (Junner & Bates 1995). 

The only major river near the Borehole Disposal Facility 
(BDF) site is River Onyasia, located about 1.3 km from the 
proposed facility and drains southwards through Achimota 
village to Accra. The Onyasia River has a depth of 0.6 m, a 
width of 6.8 m and a measured velocity of 0.8 m/s (2.5 x 
10 7 m/y). The broad valley of the Onyasia River flanks the 
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site on its eastern margin with swampy conditions generally 
found north-east of the site. Surface run-off in this area is 
very low, however, after heavy storms there is flow of water 
over the horizon below the top-soil (Darko et al. 1995). 

Seismic surveys conducted in the area mapped out 
two weak lines suspected to have been caused as a result 
of faults or fractures. However, these have not been used 
in numerical modelling to follow and predict the migra- 
tion of radionuclides at the site in case a seismic acci- 
dent occurs (Essel et al. 2011). 



Methodology 

Model definition 

The hydrogeologic setting for groundwater flow at 
steady state shall be described by the conceptual model 
shown in Figure 2. Groundwater moves from a higher 
topographic surface to a lower topographic surface. The 
water table is a free surface across which there is vertical 
recharge, denoted by R. The base of the aquifer is con- 
sidered impermeable. 



Scenario development 

In order to predict radionuclide release, the engineering 
barriers are assumed to fail in the event of a seismic activity. 
The system is thus, simplified into a two-dimensional 
conceptual model as shown in Figure 2. 

s an Open Access article distributed under the terms of the Creative Commons 
org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction 
s properly credited. 



Yeboah et al. SpringerPlus 2014, 3:155 
http://www.springerplus.eom/content/3/1/155 



Page 2 of 1 1 




The radionuclides are initially confined in the canister 
until a seismic activity occurs, leading to a crack of the 
barrier system such that the radionuclide inventory is re- 
leased into the groundwater which is the major transport 
medium. It is assumed in the calculations performed 
that radionuclides start being released from the canister 
30 years after closure of the repository. 

Numerical illustrations 

A two-dimensional numerical model was developed 
using Comsol Multiphysics (ver.3.4) similar to the pro- 
posed model in Figure 2. The lithology of the system 
was characterized by a porosity of 0.35, a transverse 
dispersivity of 0.005 m, a longitudinal dispersivity of 
0.5 m and a hydraulic conductivity ranging from 10" 15 
to 10" 5 m/s. 

The time-dependent solute transport equation as 
described by equation (1): 



° s (S) + vVC = V-pxVC)-AC 



(1) 



Here, the dispersion tensor D L defines solute spreading 
by mechanical mixing and molecular diffusion. Equa- 
tions for the tensor entries are: 



v 2 v 2 
Da = a Lr L 7 + a T T L 7 + D n 

M \ v \ 

v 2 v 2 
Djj = a L f- + a T - J - + D m 
\v\ \v\ 



(2) 



D u = Dh 



/ \ v i v i 

(a L -a T )—^ + D n 



v 2 + v 2 



Where D w Djj are the principal components of the 
dispersion tensor based on the Darcys velocity, and 
Dji are the cross terms of the dispersion tensor, the 
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Figure 2 Conceptual model. 
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subscript "L" denotes longitudinal dispersivity, "T" the 
transverse dispersivity. v is the magnitude of the Darcys 
velocity vector, D m represents the effective molecular 
diffusion in a saturated porous media and i, j are the 
spatial coordinates. 

Governing equations 

Fluid flow: assumptions 

■ On the scale simulated, the fracture system behaves 
as an equivalent porous medium 

■ The groundwater flow is assumed to be homogenous 
and subject to recharge 

■ Additionally, groundwater flows under steady state 
conditions. This means that, the velocity of flow is 
considered not to change with time since groundwater 
flow is naturally a slow process. 

Fluid flow: domain equations and boundary conditions 

Steady groundwater flow is expressed with a conservation 
equation formulated with Darcys law and expressed as: 



(3) 



Where K is hydraulic conductivity (m/s); v t is the 
Darcy velocity (m/s); x t is the spatial distance (m) and H 
is the hydraulic head (m); Hydraulic head is a function 
of pressure and gravitational potential and it is defined: 



H = h p +y 



(4) 



Where, h p is pressure head (m); and y is elevation (m). 
The equations for groundwater flow and solute transport 
are linked by: 



V(/CV/f)+£ = 0 
Where R is the recharge rate {mis). 



(5) 



Figure 3 shows the boundary conditions for the ground- 
water flow problem and stated by equations (6-8). 

A zero flux Neumann condition represents the symmetry 
boundary at x = 0 m and the impermeable boundary at 
y = 0 m as follows: 



n.(KVH) =0 



(6) 



Hydraulic head is specified at x = 14 m with a Dirichlet 
condition: 



H(x,t)=H 0 



(7) 



The annual average precipitation rate is recorded as 
800 mm/y equivalent to 2.537 x 10" 8 m/s. The water 
table receives about 10% of this value. From a topo- 
graphic height of y = 500, an inflow boundary condition 
that satisfies Neumanns condition is assumed and is 
written in the form: 



n-(KVH) =R 



(8) 



Solute transport: assumptions 

Transport of radionuclides is based on the following 
assumptions: 

■ Radioactive decay is the only reaction considered in 
the model. It is assumed to occur throughout the 
model in the liquid phase; 

■ For the purposes of this work, no gaseous release is 
considered; 

■ Transport of radionuclides is assumed to occur in the 
saturated zone; 

Solute transport: domain equations and boundary 
conditions 

For the transport equation, the initial and boundary con- 
ditions for solute transport are shown in Figure 4. In this 
work, it is also assumed that the contaminant source is 
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Figure 3 Boundary conditions for groundwater flow. 
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Figure 4 Solute transport. 



continuous over a period of thirty years. Dirichlet condi- 
tions are used at the water table, where C(x, h, t) = 0, 
except for the segment 4m < x < 8m in which: 

C = C 0 0 < t<t 0 (9a) 

C = 0 t>t 0 (9b) 

Where C 0 is the relative concentration through year 
30(£ 0 ). The Dirichlet condition at the left boundary is: 

C(0,y,t) = 0 (10) 

The Neumann condition defining the zero flow 
boundaries is: 

-n.(-6 8 D L VC)=0 (11) 

Finally, the initial condition indicating that the subsur- 
face on the site is free of contaminants at the beginning 
of the simulation is: 

C(x,y,0)=0 (12) 
Results and discussions 

Numerical simulations were considered for Cobalt-60 
(short half-life), Cesium-137 (medium half-life) and 
Americium-241 (long half-life) migrating from the re- 
pository through the subsurface environment for a 2- 
dimensional cross-section along the streamlines of 
groundwater flow. These considerations were made be- 
cause of the high amount of activity contained in the 
waste packages and the variability of their half-life as 
shown in Table 1. 

Computer simulations 

The models describe the steady-state fluid flow and fol- 
lows up with a transient solute transport simulation. 



Two partial differential equations (PDE) were solved for 
and these were assigned in separate mathematics inter- 
faces in Comsol Multiphysics (version 3.4). The first par- 
tial differential equation (PDE) is stationary and it finds 
a solution to the Darcy velocity while the second partial 
differential equation (PDE) is time-dependent and finds 
a solution to the solute transport equation. 

To test the effect of heterogeneity on solute behavior, 
a two-dimensional numerical transport model was cre- 
ated to investigate solute transport under two hydrogeo- 
logic conditions: 

(1) homogeneous hydraulic conductivity in porous 
subsurface medium and 

(2) heterogeneous hydraulic conductivity in a fractured 
medium. 

Table 2 describes the physical parameters used in the 
saturated homogenous porous medium for the numer- 
ical model of flow and transport. Except for increasing 
the hydraulic conductivity within the fracture, all other 
parameters are the same as in the homogenous model. 

Model simulation results 

Throughout the models, the amount of contaminant is 
shown by the colour bar. The activity concentration de- 
gree is indicated by the various colours, with red indicat- 
ing an intense concentration. 

Evolution of 60 Co 

Radionuclide movement (streamline plot) and concen- 
tration estimates (surface plots) of 60 Co simulated at 
various times are illustrated in Figures 5a to 6. 60 Co has 
a half-life of 5.27 years and an initial inventory of 2.22 x 
10 8 Bq/m 3 . From the simulations carried out, it was 
evident that, for a homogenous porous medium in a 
low conductivity zone, 60 Co contaminant source diffuses 
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Table 1 Ghana's inventory of disused sources (Essel 
etal. 2011) 

Disused low dose sources 



Radionuclide Total initial Application Form Unit 
activity (Bq) pieces 



Cs-137 


5.66 X10 12 


Level gauges 


Sealed 


30 


Co-60 


1 .75 x 1 0 6 


Non-Destructive 
I esting (NU I ) 


Sealed 


2 


LS- I i//LO-oU 


4.uy X 1 (J / 
4.90 x 10" 1 




Sealed 


2 


Cs-137/Am- 
241 


3.70 x 10 11 / 
1.85 X10 12 




Sealed 


3 


Cs-137/Am 
241 :Be 


3.00 X10 1 / 
1.80X10 3 


Nuclear gauges 


Sealed 


1 


Am-241 


3.50 X10 1 


Smoke detectors 


Sealed 


105 


Sr-90 


1 .25 x 1 0 4 


Thickness gauges 


Sealed 


33 


lr-192 


2.26 X10 6 


NDT 


Sealed 


1 


Cd-109 


6.66 X10 2 


Research 


Sealed 


6 


Am-241 


1.67X10 3 


Nuclear gauge 


Sealed 


1 


1-131 


6.21 X10 9 




Unsealed 


2 


Cf-252 


2.22 X10 10 




Sealed 


2 


Ra-226 


7.03 X10 9 




Sealed 


19 



needles 

H-3 3.70 x10 7 Nuclear gauge Unsealed 2750 L 

(liquid) 

C-14 2.60X10 7 - Unsealed 7empty 

(gas) cylinders 

Disused high dose sources 

1 

1 
1 
1 

2 
3 
1 
1 

2 
4 
5 

3 
4 
12 



with the groundwater which has a flow velocity of 
3.834 x 10" 8 m/s. It was observed that, the radionuclide 
migrates in the same direction as the groundwater flow 
velocity while at the same time undergoing radioactive 
decay. However, when the hydraulic conductivity was 



Table 2 Physical parameters used in the homogenous 
model for flow and transport 



Parameter 


Value 


Description 


Units 


R 


2.537e-8 


Vertical recharge 


m/s 


N 


0.15 


Effective porosity 


- 


alpha_L 


0.5 


Longitudinal dispersivity 


m 2 /s 


alpha_T 


0.005 


Transverse dispersivity 


m 2 /s 


K1 


3.17e-5 


Hydraulic conductivity 


m/s 




2.22e8 (Co-60) 






C_in 


5.66e12 (Cs-137) 


Initial activity concentration 


Bq/m 3 




3.5e7 (Am-241) 








2.78e6 (Co-60) 






D_m 


2.54e-9 (Cs-137) 


Effective diffusion coefficient 


m 2 /s 




3.17e-9 (Am-241) 








5.27 (Co-60) 








30 (Cs-137) 


Radioactive half-life 


Years 




432 (Am-241) 







increased (Figure 6) to account for heterogeneity within 
the system; it was observed that 60 Co radionuclide 
plume does not rise towards the surface. Instead, the 
plume is diluted into deeper groundwater flow systems 
as it decayed away. This may be attributed to the short 
radioactive half-life as well as the high conductivity zone 
exhibited by the fractures. 

Figure 7 shows the evolution of 60 Co contaminant 
concentration as a function of time. Here, the diffu- 
sion process takes place much longer which is ac- 
counted for by the widening of the peak hence the 
radionuclide plume concentration appears broad in 
both models. For a homogenous porous medium, the 
concentration increase steadily until it reaches a max- 
imum and then decreases exponentially. The high 
value of the initial activity concentration together with 
the diffusion process taking place, explains why the 
observed radionuclide plume extends over a very long 
range of time. 

Evolution of 137 Cs 

137 Cesium has a half-life of 30 years and an initial activ- 
ity concentration of 5.66 x 10 12 Bq/m 3 . Evolution of 
137 Cs contaminant source are shown in Figure 8a and b 
at various times (t = 250 years to t = 1000 years). In this 
case, 137 Cs is migrating in a low conductivity medium 
causing the radionuclide to be deposited onto sediments 
which sinks steeply downwards through the ground- 
water flow system. The contaminant source diffuses with 
the flowing groundwater which has a measured flow 
velocity of 3.834 x 10" 8 m/s. The flow velocity readily 



Co-60 


2.78 x10 s 


Co-60 


1.85 x10 s 


Co-60 


2.22 x10 s 


1-129* 


4.25 X10 10 


Fe-59 


2.22 X10 4 


Co-57 


1.11 X10 2 


Zn-65 


n3.70x 10 2 


Sr-89 


4.77 X10 3 


TI-204 


7.40 x 10" 1 


P-32 


1.18x10 s 


S-35 


9.25 X10" 6 / 




ml 


Ca-45 


1 .85 x 1 0 2 


Na-22 


3.7x10° 


ln-113 m 


2.22 X10 3 



Gamma cell- Sealed 
research 

Teletherapy Sealed 

Food irradiator Sealed 

sealed 

sealed 

sealed 

sealed 

sealed 

sealed 

Research Unsealed 
Unsealed 

Unsealed 
Unsealed 
Unsealed 
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moves with 137 Cs radionuclide in water causing contam- 
ination of groundwater. 

Figure 9a and b show the influence on radionuclide 
migration when the hydraulic conductivity is in- 
creased. After 40 years (Figure 9a) a fraction of its 



initial inventory would have undergone radioactive 
decay. Therefore in the presence of fractures, the in- 
crease in the hydraulic conductivity would provide 
preferential pathways for the flow velocity leading to 
an increase in mobility for the migrating radionuclide. 




Distance (m) 

Figure 6 A plot of 60 Co showing the effect of increased conductivity after 100 years. 
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The flow in this case causes the radionuclide plume to 
rise towards the surface which may finally discharge 
into wells, rivers and lakes endangering the environ- 
ment, population and biota. Figure 9b shows that, 
radioactive decay, diffusive flux and the high hydraulic 
conductivity continues to have an effect in the direc- 
tion and migration of the radionuclide source. These pa- 
rameters tend to reduce the plume concentration as well as 
cause the radionuclide source to move away from the sur- 
face towards deeper groundwater flow systems contaminat- 
ing groundwater resources. 



Figure 10 shows the evolution of Cs contaminant 
concentration as a function of time for the homogenous 
and heterogeneous model. For 137 Cs, the radionuclide 
plume concentration appears broad suggesting that, the 
diffusion process takes place much longer which is 
accounted for by the widening of the peak. Also, the con- 
centration increase steadily until it reaches a maximum 
and then decreases exponentially. The variation in the 
curves is a result of the presence of fractures which creates 
preferential pathways accounting for the high mobility of 
the radionuclide in the heterogeneous model. 




Yeboah et al. SpringerPlus 2014, 3:155 
http://www.springerplus.eom/content/3/1/155 



Page 8 of 1 1 




Evolution of Am 

The simulation was implemented for 241 Am contamin- 
ant source with an initial activity concentration of 3.5 x 
10 7 Bq/m 3 , a half-life of 432 years and a flow velocity of 
3.534 x 1(T 8 m/s. Figure 11a and b shows evolution 
of 241 Am radionuclide activity concentration at various 
times in a homogenous porous medium having a low hy- 
draulic conductivity of 3.17 x 10 5 m/s. From the 



simulations, Am was observed to sink steeply down- 
ward into the groundwater flow system after 200 years. 
At 250 years (Figure 11a), it continued to steeply decline 
into the groundwater flow system until about 1000 years 
(Figure lib). Given that the value of the flow rate is 
low due to the low value of the hydraulic conductivity, 
transport of the radionuclide was mainly dominated by 
diffusion. 
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Figure 12a to b shows the results obtained with an in- 
crease in hydraulic conductivity and its influence on the 
migration of 241 Am radionuclide. From the simulations, 
groundwater moves with a flow velocity of 13.606 m/s. 
After 60 years (Figure 12a and b), the radionuclide 
source will not have undergone any appreciable decay 
process because of the long half-life (t 1/2 = 432 years). 
Hence, majority of the radionuclide would appear to 
move towards the surface due to the high conductivity 
layer. In addition, lateral gradients can cause rapid rise 
of water-table leading to the possibility that the move- 
ment of 241 Am attributable to diffusion is faster than the 
rate at which groundwater is moving. The reduction in 
the initial activity concentration of 241 Am from 3.5 x 
10 7 Bq/m 3 to 3.256 x 10 4 Bq/m 3 can be attributed to 
mass loss by the diffusive flux. 

Figure 13 shows the evolution of 241 Am activity con- 
centration as a function of time for a homogenous and 
heterogeneous model. The nature of the peak in both 
models suggests that, advection influences the distri- 
bution moderately with diffusion being the dominant 
transport mechanism. It is noticeable that, the peak 
concentration does not reach the initial activity value 
because of mass loss by radioactive decay. Similarly, 
the concentration reaches a maximum and then de- 
creases exponentially. 241 Am activity concentration ex- 
tends over a very long range of time due to the high 
value of the initial activity, the long half-life and the 
diffusion process taking place. 

Conclusions 

The migration of three radionuclides namely, 60 Co, 
137 Cs and 241 Am have been simulated using a two- 
dimensional finite element numerical model code 
(Comsol Multiphysics). 



Neglecting heterogeneity, simulated results showed 
that, all three radionuclides ( 60 Co, 137 Cs, 241 Am) within 
the low conductivity medium sunk steeply downward 
into the groundwater flow system by diffusing into the 
flowing groundwater. This caused the flow velocity to 
move readily with the radionuclide source causing con- 
tamination of groundwater resources. 

In the presence of fractures, preferential pathways 
were created which gave rise to a rapid increase of the 
water-table and this caused the flow velocity to sweep 
the radionuclides with medium ( 137 Cs) to long ( 241 Am) 
half-life toward the surface endangering human popu- 
lation, the environment and biota. 

For 60 Co, the plume was not noticeably seen at the 
surface even in the presence of high hydraulic conduct- 
ivity but was rather diluted into deeper groundwater 
flow systems as it decayed away. This was attributed to 
the short radioactive half-life. 

The results obtained showed contamination to be 
more sensitive to variations in hydraulic conductivity 
as a result of the heterogeneity considered within the 
domain. However, impact on groundwater was still 
inevitable. 

Recommendation 

It is recommended that, proper structural geological map- 
ping including the use of stereograms should be made to 
be able to determine the fractures before radioactive waste 
is disposed of in an area. 
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